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Ph ■ Abstract 



A permanental vector is a generalization of a vector with com- 
ponents that are squares of the components of a Gaussian vector, 
in the sense that the matrix that appears in the Laplace trans- 
form of the vector of Gaussian squares is not required to be either 
symmetric or positive definite. In addition the power of the deter- 
minant in the Laplace transform of the vector of Gaussian squares, 
'sj- , which is -1/2, is allowed to be any number less than zero. 

It was not at all clear what vectors are permanental vectors. 
In this paper we characterize all permanental vectors in R^ and 
give applications to permanental vectors in i?" and to the study 
of permanental processes. 



1 Introduction 



5r , An /3-permanental vector 9 := {9i, . . . ,9n}, is an i?" valued random 

variable with Laplace transform 




E\exp{-y ^aA\ l=TrT^f^' (1-1) 
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where / is the n x n identity matrix, a is the diagonal matrix with 
ctj^j = Oj, Oj G -R+, I < i < n, and T = {r(z, j)}"^-^^ is an n x n matrix, 
/3 > and |/ + ar| > for all a G -R" . (This last requirement implies 
that detr > 0.) The fundamental question here is: For what n x n 
matrices is the right-hand side of (11 .ip a Laplace transform? 

There are very well known cases in which right-hand side of (II. ip is 
a Laplace transform. When /3 = 1/2 and F is symmetric and positive 
definite, 

9=iGl/2,...,Gl/2) (1.2) 

where {Gi, . . . ,Gn) is a Gaussian random variable with covariance F. 
(We sometimes refer to a vector like ^ as a vector of Gaussian squares.) 

The innovation in the question posed here is that we consider all 
P > and do not require that F is symmetric or positive definite. 

Before we attempt to answer this question it is important to note that 
the matrix F is not unique. If D is any diagonal matrix with non-zero 
entries we have 

\I + aT\ = \I + aDTD-^l = \I + aDT'^D'^\, (1.3) 

for all diagonal matrices a. The matrix F is said to be diagonally equiv- 
alent to F' if F' = DTD^^ for some diagonal matrix D with non-zero 
entries. For a very large class of irreducible matrices F, it is known 
that the class of diagonally equivalent matrices are the only sources of 
non- uniqueness; see [7]. 

Sometimes one can take D to have diagonal entries ±1. Such matrices 
are called signature matrices. It is obvious that if S* is a signature matrix 
then S = S^^. We also note that (11.31) may hold with D = I, the identity 
matrix even when F 7^ F'. For example, if F and F' are n x n matrices 
with the same diagonal elements and all zeros below the diagonal, then 
(11.30 holds with D = I. In this case we say that F and F' are effectively 
equivalent. (We also note that we sometimes refer to F as a kernel for 
6.) 

R. Bapat and R. C. Griffiths [HH], (see, also [H Chapter 13]), com- 
pletely describe the vectors of Gaussian squares for which, (ll.ip is a 
Laplace transform for all /3 > 0. They do this in solving a classical prob- 
lem posed by P. Levy: When is a vector of Gaussian squares infinitely 
divisible? The answer is: 



Proposition 1.1 A vector of Gaussian squares is infinitely divisible if 
and only if the covariance matrix T is diagonally equivalent to an M- 
matrix. 

A matrix A = {ajj}, is said to be an M matrix if 

(1) Ojj < for all i ^ j. 

(2) A is nonsingular and A^^ > 0. 

Strictly speaking knowing that a vector of Gaussian squares is in- 
finitely divisible only asserts the existence of the Laplace transform for 
(3 = k/{2n), for all integers k,n > 1. However the proof of Proposition 
O shows that dH]) holds for all /3 > 0. 

There are permanental vectors with kernels that are not diagonally 
equivalent to symmetric matrices. Eisenbaum and Kaspi, P, Lemma 
4.2] recognize that the Bapat-Griffiths sufficient condition for infinite 
divisibility in the case of symmetric kernels also works for non-symmetric 
kernels. 

It is well known that positive definite symmetric matrices may or may 
not have an inverse that is diagonally equivalent to an M-matrix. On 
the other hand when /3 = 1/2, (11. ip is the Laplace transform of a vector 
of Gaussian squares. Based on these observations we divide the class of 
kernels F of permanental vectors into three categories. 

1. r is diagonally equivalent to a symmetric positive definite matrix. 

2. r~^ is diagonally equivalent to an M matrix. 

3. r is not in class 1. or class 2. 

Note that we use the expression positive definite to include what is 
sometimes called positive semi-definite. Also, we emphasize that classes 
1. and 2. are not disjoint. 

There is ample reason to think that there is an abundance of examples 
of kernels in class 3. One should be able to take a symmetric positive 
definite matrix not in class 2. and alter its off diagonal elements very 
slightly. One then might expect that the altered matrix would be in 
class 3. We worked for a long time to find an example of a kernel of 



a permanental vector in class 3. but were not successful. We then set 
out to fully characterize 3x3 matrices that are kernels of permanental 
vectors. The main result of this paper is the following theorem which 
states that for permanental vectors in R^, class 3. is empty: 

Theorem 1.1 A 3 x 3 matrix that is the kernel of a permanental vector 
in R\, belongs to class 1. or class 2., or both. 

This result also applies to permanental vectors in i?" in the sense 
that if 6* = {^1, . . . , On} is a permanental vector in i?" then any three 
components of ^ is a permanental vector in R\. 

Another consequence of Theorem 11.11 is that if a kernel of a perma- 
nental vector in R\ is not diagonally equivalent to a kernel with positive 
entries then it is the kernel of a vector of Gaussian squares. 

A permanental process {P(t),t G T} is a stochastic process with fi- 
nite dimensional distributions that are permanental vectors. Eisenbaum 
and Kaspi study permanental processes in [3]. Roughly speaking they 
show that the the potential density of a Markov process is the kernel of 
a permanental process. (When this is the case we say that the perma- 
nental process is associated with the Markov process.) In fact they show 
that permanental processes are the missing link that allows the Dynkin 
Isomorphism Theorem to be extended to the local times of Markov pro- 
cesses that are not symmetric, [21 Corollary 3.5]. There are several in- 
timate connections between permanental processes with a kernel that is 
the potential density of a Markov process and the Markov process itself. 
In P the permanental process is shown to be the loop soup local time of 
the Markov process. 

Permanental processes are introduced by Vere- Jones in [10]. In [TOl 
Proposition 4.5] he gives necessary and sufficient conditions for (11. ip to be 
the Laplace transform of the vector (^i, . . . , ^„) in terms of the modified 
resolvent matrix 

r, :=r(/ + rr)-^ (i.i) 

where r > 0, and F is the matrix in (II. ip . 

Proposition 1.2 nO|. Proposition 4-5] For U.l\) to represent the Laplace 
transform of a non-negative random vector it is necessary and sufficient 
that for all r > 



(i) Tr exists and is (5-positive definite, 
(ii) det(/ + rr) > 0. 

Furthermore, given F^, Proposition II ■2[ (i) may hold for some values 
of /3 but not for others 

(Item (ii) is equivalent to: All the real, non-zero, eigenvalues of F are 
positive.) 

There is no point in giving the very complicated definition of /3- 
positive definite here. One can refer to [10] or to [31 |6] where it is 
repeated. It seems almost impossible to verify Proposition 11.21 (i) un- 
less all the entries of the matrix F^ are greater than or equal to zero, 
in which case (i) holds for all /3 > 0. In [Sj Theorem 3.1] Eisenbaum 
and Kaspi point out that this is is the case when F is associated with a 
Markov process and that Proposition [L2] (ii)-, also holds for these kernels. 

If all the entries of the matrix F^ are not greater than or equal to zero 
verifying Proposition 11.21 (i) necessitates examining an infinite sequence 
of increasingly larger matrices derived from F^. Otherwise we know no 
sufficient condition for the existence of a permanental vector that might 
be in class 3. One is given in [TOj Proposition 4.6], but it is not correct. 
We discuss this in Remark 18.11 

There is a potentially important application of Proposition 11.21 if one 
can figure out how the verify (i). If the kernel of a permanental vector 
is in class 2., (11. ip is a Laplace transform for all /3 > 0. If the kernel 
of a permanental vector is in class 1. and not in class 2. then we only 
know that (II. ip is a Laplace transform for /3 = 1/2, and trivially, for all 
/3 = A;/2, for integers k > 1. Possibly there exist other values of /3 > 
for which (II. ip is a Laplace transform. Applying Proposition II. 2| which 
depends on /3, would answer this question. 

There are many other interesting applications of Theorem 11.11 The 
next result answers a question that started our interest in 3-dimensional 
permanental vectors. We point out in the beginning of this Introduction 
that the univariate marginals of a 1/2-permanental process are squares 
of normal random variables. It also follows from (II. ip that pairs (6*^, 9j), 
of a 1/2-permanental process, are equal in law to (6*1/2, 6*1/2), where 



(Gj, Gj) is a Gaussian vector with covariance matrix 
f - 






:i.5) 



(See [9l Lemma 3.1]). It follows from this that 



T?m\ ^(^'^) A fa a\ r(i,j)r(j,i) 

^(fe'i) = and cov|6'i,0j| = . (1.6) 

Therefore, if 

r(z,j)r(j,z)=0 Vl<z^j<n (1.7) 

the components of a 1/2-permanental process are pairwise independent. 
Actually, (I1.7P is a necessary and sufficient condition for the compo- 
nents of any /3-permanental process to be pairwise independent. This 
is because in this case the determinant of F is a product of its diagonal 
elements and the right-hand side of (II. ip 

|/ + aTf = \I + a,T{i,t)f\I + a,T{j,j)f. (1.8) 

We also know from [THI, bottom of page 135] that for any /3 permanental 
process 

cov{e,,e,} = PTit,j)Tij,t). (1.9) 

If r is symmetric and positive definite and 9 = {Gl/2, ... ,01^/2), 
where {Gi, . . . ,Gn) is a Gaussian random variable with covariance F, 
with T{i,j) = F(j, i) = 0, then the components of 9 are independent. We 
asked ourselves the following question: "For a general /3-permanental 
vector 9, that is not the square of a Gaussian vector, does (I1.7P imply 
that the components of 9 are independent?" The answer is yes. We 
prove: 

Theorem 1.2 Let 9 be an n- dimensional f3 -permanental vector with pair- 
wise independent components. Then the components of 9 are indepen- 
dent. 

It is clear that when /3 = 1/2 Theorem 11.21 implies that the only 
permanental vectors with independent components are those with com- 
ponents that are squares of independent Gaussian random variables. 



The next result deals with a function that appears in sufficient con- 
ditions for the continuity of permanental processes in ^, 

d{x, y) = {t{x, x) + T{y, y) - 2 (r(x, y)T{y, x))''^) ^'^ . (1.10) 

If r(a;, y) = T{y, x) is the covariance of the Gaussian vector {G(t),t G T}, 
then 

dix,y)={EiGix)-Giy))Y\ (l-H) 

which is a metric on T. However, if r(a;, y) is the kernel of a permanental 
process and r(a;, y) ^ T{y, x) it was not clear whether or not d{x, y) is a 
metric on T. We can now say that even in this case d{x, y) is a metric 
on T. 

Corollary 1.1 Let {P{t),t E T} be a permanental process with kernel 
r{x,y). The function {d{x,y),x,y G T} in Iil.l0\) is a metric onT. 

In Section [2] we give many properties that are necessary for an n x 
n matrix to be the kernel of a permanental vector. In Section [3] we 
obtain an interesting property of the eigenvalues of 3 x 3 positive definite 
symmetric matrices that plays a critical role in the proof of Theorem ll.il 
The proof of Theorem 11.11 uses completely different methods when the 
off diagonal elements of the kernel are all negative or all positive. These 
cases are considered separately in Sections S] and O SectionsE]and[7]give, 
repectively, the proofs of Theorem 1 1 . 2 1 and Corollary ll.il Finally, because 
the fact that kernels in class 2. are kernels of permanental vectors is so 
important in this paper, we give an outline the proof, essentially showing 
what changes are necessary in the proof in [8], Theorem ] , which is given 
for symmetric kernels. This result is given in [31 Lemma 4.2]. The proof 
involves probabilistic considerations. Since Theorem 11.11 is only for 3x3 
matrices it seems appropriate to give a proof for finite matrices involving 
only linear algebra. 

We are grateful to Professor Jay Rosen for many helpful comments 
and discussions 



2 Preliminaries 

If ^ is a permanental vector in i?" then any subset of its components, say 
of p components, is a permanental vector in i?!L. For p = 2, the Laplace 



transform of the vector {6i, 6j} takes the form 

^(exp(^-iM, + a,0,)^^ (2.1) 

: {l + aiT{i,i)) + ajT{j,j) 



\i + ar\ 

+aiaj{T{i,i)T{j,j) -T{iJ)T{j,i))y'^ . 
Taking ctj = aj sufficiently large, this implies that 

r{z,z)r{j,j)-r{t,j)r{j,t)>o. (2.2) 

If we also set aj = in (11 .ip we see that for any z G n 

T{i,i)>0. (2.3) 

In addition, by [10, Proposition 3.8], for any pair i,j G T 

T{t,j)T{j,t)>0. (2.4) 

In the next lemma we show that there are many transformations of 
kernels of permanental processes that give other kernels of permanental 
processes. 

Lemma 2.1 Let A be a kernel of a P -permanental vector 9 = {6i, . . . , 
9n)- Let Ui and U2 he diagonal matrices with non-zero diagonal entries 
uf , i = 1, . . . ,n, j = 1,2, for Ui and U2 respectively, with the property 
that u\ u\ > 0, i = l,...,n. Then U1AU2 is the kernel of the /3- 
permanental vector {u[ u[ 61, . . . ,Un Un On)- 

Proof Since 9 is an i?" valued random variable so is {u\ u[ 9i, . . . , 



Un Un On)- The Laplace transform of {u\ ui 61, . . . ,Un Un On) is 
E [ exp ( - J2 cy^{u^^uf^9A ] = E ( exp ( - ^^(a.nf )«, 



4\^^^^.)11 = i^fexpf-g^--^^^-^^^^ 



\I + {aU2Ui)A\~^ (2.5) 

\U2{I + a{UiAU2)U2T'' 
\I + a{UiAU2)\-^ - 



n 



Example 2.1 We note two cases. Let f/ be a strictly positive diagonal 
matrix 

1. When 7 + 7'=0 

Elexpl-J2 "*^i ) ) = 1^ + aU^^U-^\~^- (2-6) 

2. When 7 + 7'=! 

E (exp ( -^a^Mi^i j I = \I + aU^AU^^-^^\-^. (2.7) 

In particular, 

E [exp j-^a,M,^,] I =\I + aU'^^AU^/^\-^. (2.8) 

Remark 2.1 It is easy to see that for a',b',c' strictly positive, the two 

matrices 

' 1 a' c' \ / 1 -a' -d \ 

a' I V \ and -a' 1 6' . (2.9) 

d V I j \ -d V 1 / 

are diagonally equivalent to each other. Similarly 

I -a' d \ /I -a' -d \ 

-a' 1 y and -a' 1 -V (2.10) 

d y I j \ -d -V 1 / 

are diagonally equivalent to each other. It should also be clear that these 
observations hold if any two of the three pairs of entries in (12.91) are taken 
to be negative and if any pair of entries in (12.101) is taken to be negative. 
Because of the observations in the previous paragraph, when we con- 
sider whether a 3 x 3 matrix is the kernel of a permanental vector we 
need only consider those matrices with all positive off-diagonal elements 
or all negative off-diagonal elements. (We consider to be both positive 
and negative.) 



Consider the matrix 




E= \ a2 1 h \ . (2.11) 



The next elementary lemma is very useful. We leave the proof to the 
reader. 

Lemma 2.2 // the off diagonal elements of the matrix E in Ii2.11\) are 
either all strictly positive or all strictly negative then E is diagonally 
equivalent to 

(1 a c \ 
a 1 b[ \ . (2.12) 

c b', 1 / 

where a^ = 0102, c^ = C1C2 and h'ih'2 = &1&2 (I'fT'd the signs of a, c, b'l, h'2 are 
the same as the signs o/ai, 02, Ci, C2, &i, &2- 

By (12. 4p the kernel F of a permanental vector has the property that 
r(i, j)r(j, i) > 0. Therefore, if E is the kernel of a permanental vector 
in i?i^, ai, a2 are either both positive or both negative, and similarly for 
61,62 and ci,C2. 

Lemma 2.3 Suppose that the matrix E is the kernel of a permanental 
vector in R\. Then if 

aibiCi = a2b2C2 (2-13) 

it is diagonally equivalent to the kernel 



£ = \ ±a 1 ±b \ (2.14) 




where a = (0102)"^'^, b = (6162)^ andc = (0102)^'^ and in which {^}i,2 = 
a if ai is positive and {^}i,2 = —a if ^i is negative, and similarly with 
respect to 61 and ci . 

In particular this lemma holds when both sides of ^2.1^) are equal to 
zero. 



10 



Obviously, £ is the covariance of a Gaussian vector. 
Proof When fl233|) holds and aihiCi ^ 

aibiCi + a2b2C2 = 2(aia2&i&2CiC2)^''^ (2-15) 

It is easy to see that E and £ are diagonally equivalent. It is also easy 
to see that E and £ are diagonally equivalent if, say, 61 = 62 = 0. 

Finally, it is also easy to see that if ai = 62 = 0, ii^ is effectively 
equivalent to 

/ 1 (CiC2)l/2\ 

£' = \ 1 . (2.16) 

V (ciC2)^/2 1 / 

D 

We also use the following lemma which is [Gj Lemma 4.5]. 

Lemma 2.4 Let 

(u a c \ / u ai C2 \ 

a V b \ , B= I a2 V bi \ . (2.17) 

c b w J \ ci b2 w J 

where aia2 = a^ , 6162 = b'^ ,ciC2 = c^. Suppose that ^ > 0. If B~^ is an 
M matrix then A^'^ is an M matrix. 

Remark 2.2 By definition an M-matrix is invertible. Therefore |v4| > 
0. We also note that if B~^ is an M matrix then for any diagonal matrix 
D with strictly positive entries, {DBD~^)~^ is an M-matrix. 

The next observation is used often in this paper. 

Lemma 2.5 Let $(ai, . . . , a„) be the Laplace transform of an i?" valued 
random variable. For any 1 < k < n set aj = Uj, where Uj > 0, 
k < j < n. Then 

t,„..,K,...,«>)- '''';^--°"'"'^'--''"' (2.18) 

$(0, ...,0,1ifc+l,...,M„) 

is the Laplace transform of an R^ valued random variable. 

Furthermore, if^{ai, . . . , a„) is the Laplace transform of an n-dimen- 
sional permanental vector, $(„^fc)(ai, . . . , ak) is the Laplace transform of 
a k-dimensional permanental vector. 

11 



Proof Since $(ai, . . . , a„) is a completely monotone function on i?" it 
follows that ^(^n,k){<^i, ■ ■ ■ , ctfc) is a completely monotone function on R^, 
satisfying $(„^fc)(0, . . . , 0) = 1. Therefore, it is the Laplace transform of 
an R^ valued random variable. 

(This is very well known when k = 1. Lacking a suitable reference for 
general k, we note that it follows from the Extended Continuity Theorem 
for probability measures on i?^, [SI Theorem 5.22], and the argument 
in the proof of [H Lemma 13.2.2], applied to $(„_a:)(«i, • • • ,«fc), not its 
logarithm. ) 

Now suppose that $(ai, . . . , «„) is the Laplace transform of an n- 
dimensional permanental vector. This implies that 

$(«!,...,«„) = -—^-^ (2.19) 

for some n x n matrix F, and diagonal matrix a as in (II. ip . We first 
prove the second statement in the lemma for k = n — 1. Consider 

$(ai,...,a„_i,u„) (2.20) 

and the corresponding matrix / + aF, where a = (ai, . . . , a„_i, m„). 
We now show that 

|/ + 5F| = (1 + Un)\I + a(""i)F("-^)| (2.21) 

where a'-"'"^^ is the (ra — 1) x (n — 1) diagonal matrix with diagonal entries 
(ai, . . . , a„_i) and F'^""-'^^ is an (n — 1) x (n — 1) matrix with entries that 
are functions of the entries of F and «„. Since 

$(0,...,0,...,0,^„)= ^ (2.22) 

|1 + Un\P 

The equality in (I2.2ip gives (I2.18P when k = n — 1. 

To obtain (I2.2ip we note the matrix I + 5F has the same determinant 
as the matrix obtained from it by subtracting F(n, j) ^";^ . times the n-th 
column from the j-th column, for each 1 < j < n — 1. Call this matrix 
S. Note that S{n,j) = 0, j = 1, . . . ,n — 1 and S{n, n) = 1 + Un- Let S' 
denote the martix obtained by dividing the last row of 5* by 1 + m„. We 
have 

|/ + 5F| = (1 + M„)|5'|. (2.23) 

12 



To be more specific the entries of S' are 

S\z,j) = <5,, + a.fr(.,j)- ''"^^^'!^^^'''^'^ V l<^,J<n-l■ 

S'{n,j) = 0, l<J<n-l; 

S'{n,n) = 1. (2.24) 

It is obvious that we can write 

|5'| = |/ + a("-')r("-^)| (2.25) 

where T^'^"^^ is the matrix with components 

u,.nM,lY \ 1 <,,<„_!, (2,26) 

1 + Un J 

and a^^~^^ = {ai,a2, ...a„_i). We now have 
Repeating the argument above we can show that 

^ / N ${n,n-l)(ai,---,an-2,Mn-l) .^ ^„. 

$(n,n-2) ai,--.,an-2) = — T 77^ ^^ ^ 2.28 

$(n,n-l)(0, ...,0,M„_i) 

_ $(ftl,...,ftn-2,'»n-l,'»n) 
$(0, ...,0,U„_i,U„) 

since 

^ . X $(ai,...,an-2,M«-l,M«) .„ „„s 

$(n,n-l) (ai, • • • , an-2, Mn-l) = ^^77^ ' 7 (2.29) 

$(0, ...,0,M„) 

and 

^"'^ '^^ ' ' ' " 1^ $(0,...,0,u„) ^ ' 

Thus we get (12181) for A; = n - 2. 

Continuing in this way we get (12.181) for arbitrary 1 < k < n — 1. n 

We use the following necessary condition in the proof of Theorem 1 1.1[ 
It is a direct consequence of [TOl Proposition 3.8 and Proposition 4.5 with 
o"=0]. We provide a direct proof for the convenience of the reader. 

13 



Lemma 2.6 Let A = {Aij}^-^, be an n x n matrix. If A is a kernel 
of a P -permanental vector then A and all matrices obtained from A by 
multiplying its rows by non-negative numbers have a positive eigenvalue 
of maximum modulus. 

Proof Let 6 = {6i, . . . , On) be a perinanental vector with kernel A. The 
matrices obtained by multiplying the rows of A by non-negative numbers 
have the form UA where [/ is a diagonal matrix with non-negative entries 
ui . . . ,Un- Note that 

E [ exp [ - ^ aiUiOi j j = |/ + aUAl-l^. (2.31) 

Let z be a complex number and set 

f{z) = E lexp IzY,ua]] (2.32) 

" ^ kl 

k=0 

By (EUD 

n 

f{z) = \I- zUAf = n (1 - ^^p)'^ (2-33) 

p=i 

where Ap, 1 < p < n, are the eigenvalues of UA. 
Since 

we see that 

EieX) = iAl)TiP + k). (2.35) 

It follows that 

Eiy^Uiei] < maxu^n'' maxE{e^) (2.36) 

\ ■^— ^ / l<i<n l<i<n 



i^n'' max '^^ 

l<i<n l<i<n 



max ufn" max AlJ:{l3 + k). 
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Consequently, there exists a number b, such that 

k\ 



(2.37) 



This iniphes that the series in fl2.32p has a positive radius of convergence, 
which we denote by R. By fl^:^ and (105]) 



f{z) = J2^kz'' = llil-z\ 



1-/3 



(2.38) 



k=0 



p=l 



for |z| < i?. Let f > and note that lim^^ji f{v) = oo as the sum of a 
series with positive terms. 

Since the terms Bk are positive, when v = \z\ we have 



f{v) = J2^k\z\'> 



k=0 



fc=0 



kZ 



(2.39) 



This shows that if f{z) has a singularity at Zq, then /(fo) = oo, for 

^0 = \zo\- 

By fl^:^ this implies that 



max Af 

l<p<n 



1 



(2.40) 



□ 



3 Eigenvalues of 3 x 3 positive definite sym- 
metric matrices 



Lemma 3.1 Let H he the real symmetric matrix 

lac 

H= \ a I b 

c b 1 

and let p be a real diagonal matrix with entries (pi, p2, ps), where 

b c a 

Pi = -, , P2 = 7, and p3 



(3.1) 



b — ac 



c — ab 
15 



a — be 



(3.2) 



Assume none of the denominators in li3.S\) are zero. Then 

\pH - XI\ = (A - l)\p,p,p,\H\ - A). (3.3) 

In particular, X = 1 is an eigenvalue of pH of multiplicity 2. 
Proof 



\pH-I\ 



ab be 



b—ac 
ac 


b—ae 

ab 
e—ab 

ab 
a~be 


b—ae 
be 


c—ab 
ac 


e—ab 
be 


a— be 


a— be 



(3.4) 



It is easy to see that the second and third row of this determinant are 
equal to a (different) muhiple of the first row. This shows that A = 1 is 
an eigenvalue of pH of multiplicity 2. 

Since the product of the eigenvalues must equal \pH\ we get (l3.3p .D 



Remark 3.1 When detif > 0, or equivalently, when H is positive def- 
inite, H is the covariance of a Gaussian vector, say, (^1,^2,^2)- When 
p > 0, it follows from Lemma [2711 that pH is diagonally equivalent to the 
covariance matrix of the Gaussian vector 

(pl/^Gi,pfG2,Gfe3). (3.5) 

4 3x3 matrices with negative off diagonal 
elements 

In this section we prove Theorem 1 1.1 1 for matrices with negative off diag- 
onal elements. 

In the next lemma we consider the eigenvalues of a 3 x 3 matrix with 
negative off diagonal elements. 

Lemma 4.1 Let A_ he the matrix 



1 


-< 


-4 


-4 


1 


-&; 


< 


-V, 


1 



A_ = -4 1 -6; . (4.1) 
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with a'^a^i h'^h^., d^d^ all greater than or equal to and less than or equal to 
1. Assume that del A > and AijAj^i < 1. Then if A^ is not diagonally 
equivalent to a symmetric matrix there exists a diagonal matrix $, with 
strictly positive entries, such that $y4_ has only one real eigenvalue. 

Furthermore, the real part of the complex eigenvalues of $y4_ is 
greater than the real eigenvalue. 

Proof We first consider the case in which a'^, a'g, &'^, 62, c'^, C2 are all 
strictly positive, Let a'ia'2 = a^, b[b'2 = 6^ and c[d2 = c^ and let $ be a 
real diagonal matrix with entries (0i, 02, ^s), where 

and 03 = ^— . (4.2) 



b + ac c + ab a + be 

Let K be the matrix 



K= \ -a 1 -b \ . (4.3) 




We show in Lemma [3. II that 
\^K-XI\ (4.4) 

= -A^' + A^ (01 + 02 + 03) 

-A (0102(1 - a^) + 0203(1 - &') + 0103(1 - C')) + 010203|i^| 

= (A-l)2(0i0203|ir|-A). 

Since (0i + 02 + 03)/3 < 1, the second derivative of \^K — XI\ is negative 
at A = 1. Therefore, \^K — A/| has a local maximum of at A = 1. 
Consider 

H{\):=\^A_-XI\. (4.5) 

By Lemma 12.21 it suffices to take A_ to be 



(4.6) 





/ i 


—a 


— c 




-a 


1 


-bi 




V- 


-62 


1 


where 6162 = b"^- 








We have 









|ir| = 1 - a^ - 52 _ c^ - 2abc (4.7) 
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and 

\A_\ = l-a'^ -b"^ -c'^-ac{bi + b2). (4.8) 

Therefore, since 6162 = b'^ we see from the first equahty in (14. 4p that 

\^K - XI\ - |$A_ - A/| = 010203 (ac(6i + 62 - 26)) . (4.9) 

Unless, 61 = 62, the right-hand side of (14. 9 p is strictly positive. Since 
\^K — A/ 1 has a local maximum of at A = 1, |$A_ — A/| has a local 
maximum that is strictly negative at A = 1. This implies that when 
61 7^ 62; \^A_ — XI\ has only one real root. This is equivalent to the 
statement of this lemma. 

We now consider the cases in which some of the entries of A^ in (14. ip 
are zero. Considering Lemma 12.21 we can restrict our attention to the 
following matrices 

1 -a \ / 1 -a' 

-a 1 -b' \ 1 -b' 

-c' 1 / \ -c' 1 

(4.10) 
in which a, b, c, a',b', d are all strictly positive . Label them, respectively, 
J^i, 72 and J's. 

We first show that Tx has only one real eigenvalue. Similar to (14. 9 p 
we have 

\^K-\l\-\^Ty-\l\ (4.11) 

= -0203^^ (01 - A) - 010203 (2a6c - acV) 
= 0203&^A — 010203 {b'^ + 2abc — acb') . 

We assume b' is fixed. One can choose < 6 < 1 such that 6^ + 2abc — 
acb' = 0. That 6 > is elementary. That 6 < 1 follows from the fact 
that det J-'i > implies that acb' < 1. With this choice of b 

|$ir-A/| - I^J'i-A/l =02036^A. (4.12) 

Considering the graph of \^K — A/|; (see the last line of (I4.4p ). and the 
fact that the right-hand side of (I4.12p is strictly positive for A > we see 
that |$J^i — A/| has only one real root, or equivalently, that J-'i has only 
one real eigenvalue. 
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Similar arguments show that T^ has only one real eigenvalue. With 
regard to T2 we have 

\^K-\l\-\^J=2-\l\ (4.13) 



= -0203^^ (01 - A) - 0i03C^ (02 - A) - 010203 (2a&c - ad^J) 

= (0203&^ + 0103C^)A - 010203 (&^ + C^ + 2a6c - flC 6') . 

For simplicity we can take b = c. We can can choose 1 > 6 > such that 
26^ + 2ab^ - ac'b' = 0. With this choice of b 

|$K - A/| - |$J-2 - A/| = (0203 + 0103)&'A. (4.14) 

So we can use the same argument we just used to show that J^2 has only 
one real eigenvalue. 
For J-'a we consider 

I^Js-A/l-I^J'a-A/l (4.15) 

= 0i02a^A — 010203 (a^ + ab'c' — a'b'c') . 

We can choose < a < 1 so that a^ + 2ab'c' — a'b'c' = and get 

I^J's - A/| - I^J's - A/| = 0i02a^A. (4.16) 

Considering the graph of |$J^2 — A/| we see that J-3 has only one real 
eigenvalue. (Actually, in this case, it is easy to see that the three eigen- 
values of J-3 are {{abcY^^, (a6c)^/3^i±iv^| ) 

This completes the proof of the assertions in the first paragraph of 
this lemma. 

We now consider the assertions in the second paragraph of this lemma. 
Consider H{X) Suppose that 5 is a real root of this polynomial. Dividing 
if (A) by A — 5 we see that the real part of the complex roots of H{X) is 
greater than 6 if and only if 

01 + 02 + 03 > 35. (4.17) 

Suppose that a[,a2,b[,b'2,c[,c'2 are all strictly positive. It follows 
from (14.91) and the fact that 010203 det A is a root of \^K — A/|, that 
6 < 010203 det ii'. Furthermore, since \K\ < 1, to prove fl4.17p it suffices 
to show that 

01 + 02 + 03 > 3010203. (4.18) 
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Since 01 + 02 + 03 > 3(010203)^''^ and 0i, 02, 03 are all less than 1 we see 
that f l4.18p is satisfied. Thus, when a\, a!^^ b[, 62, c[, Cg are all strictly pos- 
itive and $74_ has only one real eigenvalue, the real part of the complex 
eigenvalues is greater than the real eigenvalue. It is easy to see that this 
argument also works when y4_ has the form of J^i, T2 or J-3. n 

Proof of Theorem 11.11 for matrices with negative off diagonal 
entries This follows immediately from Lemmas 14.11 and 12.61 n 



5 3x3 matrices with positive off diagonal 
entries 

Consider the matrix 



A=\ a^ 1 61 , (5.1) 




in which all the entries are greater than or equal to zero. When det A > 0, 

. / 1 ~" ^1^2 £2^2 - ai fli&i - C2 
A~^ = — bici - 02 1 - C1C2 a2C2 - 61 I . (5.2) 

' ' \ a2&2 - ci aiCi - 62 1 - 0102 

Proof of Theorem 11.11 for matrices with positive off diagonal 
entries This follows from the next lemma: 

Lemma 5.1 Let A, in Ii5. 1\) . be the kernel of a permanental vector 9 = 
{9i, 62,93). Then A is either diagonally equivalent to a symmetric positive 
definite matrix or A~^ is an M matrix. 

Furthermore, if one of the off diagonal terms of A is equal to 0, then 
A is diagonally equivalent to a symmetric matrix. 

Proof As in (I2.19p . but with F replaced by A, we have 

$(«!, 02, as) = \j ^^Mp - (5-3) 
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Therefore, as in (12.271) 

,. .. (^. ^„^ - 

|/ + a(2)r(2)|/3' 

where, by (Km 



•^(1,2) («i, «2) = I, , ^(,^r^(,^^, ■ (5-4) 



\^ a2 - t-^iCi 1 - vbib2 J 
and where f = . It follows from Lemma [2.51 that r*^^) jg ^]-^g kernel 

1 + ^3 

of a permanental vector. Therefore, by (12. 4p 

{ai - vc2b2){a2- vbiCi) >0. (5.6) 

Suppose none of the off diagonal entries of A are equal to 0. The inequal- 
ity in (15. 6 p holds for all v G (0, 1). Therefore, either 

ai > C2&2 and 02 > &iCi (5.7) 

or there exists a f G (0, 1) such that 

oi - VoC2b2 = 02 - VoCibi = 0. (5.8) 

It follows from (15.80 that ai = t'oC2&2 and 02 = foCi^i, or, equivalently, 
that 

aibiCi = a2b2C2. (5.9) 

If (15. 9p holds, it follows from Lemma [2.31 that the matrix A is diagonally 
equivalent to a symmetric matrix. 

We repeat this argument twice considering $(3,2) («i, ^3) and $(3,2) (0^2 , 
as). If (15. 8p holds then we get comparable equalities when we consider 
$(3,2)(«i, ct3) and $(3,2)(ct2, a3)- However if (15.70 holds we also get 

61 > a2C2, ^2 > fliCi, ci > a2&2 and C2 > ai6i. (5.10) 

It follows from (15. 7p and (15.100 that A is an M-matrix if it is invertible, 
or equivalently, \A\ > 0. 

However, if |y4| = 0, A does not have an inverse and the consideration 
of whether A~^ is an M-matrix is meaningless. Therefore we must show 
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that when f l5.7p and f IS.lOp hold and A is not diagonally equivalent to a 
symmetric positive definite matrix then \A\ > 0. 

We need only consider the case in which A is not symmetric. Without 
loss of generality we can consider that 




A= \ a I b \ . (5.11) 



where a,b,c > 0, a^ = aia2, b"^ = foi&25 C1C2 = c^ and Ci 7^ C2. Let d be 
such that Ci + C2 = dc. Obviously d > 2. 
Since \A\ = we have 

1 - (a^ + 6^ + c^) + abed = (5.12) 

We consider d in (15.121) as a function of a, b, c, i.e., 

d[a,b,c) = ; . (5.13) 

abc 

Note that the gradient of d{a, b, c) 

Vd{a, b, c) = — — (6c(l + a^-b^- c"), (5.14) 

ac(l + 6^ - a^ - c^), a6(l + c^ - a^ - b'^)) . 

The inequalities in (15. 7p and (I5.10p hold when A is not diagonally equiv- 
alent to a symmetric matrix. (The argument we give does not require 
that |y4| > 0.) When they hold we see that the components of Vd{a, b, c) 
are all greater than or equal to 0. For example, a > be implies that 

{1 + a^-b^- c^) > (1 + {bcf -b^- c^) = (1 - b^){l - c^) > 0. (5.15) 

Note that d{l, 1, 1) = 2. Therefore, since a, b, c are all less than or equal to 
1, d{a,b,c) < 2. This contradiction shows that here are no permanental 
vectors with |yl| = other than those with kernels that are diagonally 
equivalent to a symmetric matrix. 

To show that if one of the off diagonal terms of A is equal to 0, 
then A is diagonally equivalent to a symmetric matrix we consider (15. 6p . 
Suppose ai = then either 02 = or one of 62, C2 is equal to 0. In 
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these cases it follows from Lemma 12.31 that A is effectively equivalent 
to a symmetric matrix. Using <l>(2,3) and $(1,3) we come to the same 
conclusion for all the other ways one or more of the off diagonal terms of 
A can be equal to 0. n 

Example 5.1 It seems worthwhile to point out that there are many 
symmetric matrices with positive entries that have determinant 0. All 
3x3 symmetric matrices of the form of D in fl3.ip . with, \a\ < 1, |6| < 1 
and |c| < 1 and with |D| = have the form 

(1 sinx cosy \ 

sinx 1 sm{x±y) , (5.16) 

cosy sm{x±y) 1 J 

for any x and y which sin x, cos y, sin(x + y) or sin x, cos y, sin(x — y) are 
greater than or equal to zero. To get this we note that |D| = implies 
that 

c = ab±{{l-a^){l-b^)Y^\ (5.17) 

If we take a = sinx, and b = cosy and solve for c we get (I5.16p . 

We now ask for what values of x and y is the adjugate, (also called the 
adjoint) of S±{x,y) a singular M-matrix. (I.e. even though the matrix 
is not invertible, the adjugate has negative, including 0, off diagonal 
elements.) Referring to D, and noting (15.21) . this is equivalent to asking 
for what values of x and y are 

c> ab, a > be, and b > ac. (5.18) 

To achieve the first inequality in (I5.18P we must use the plus sign in 
(I5.17P which gives c = sin(x + y). This implies that a, which satisfies an 
analogue of (I5.17p . satisfies 

a = be - {{I - b^){l - c'')y^\ (5.19) 

To get the second inequality in (15.181) we can take c = 1 and a = b, which 
can be achieved by taking x G [0, 7r/2] and y = (vr/2) — x, or a = c and 
6=1 which can be achieved by taking x G [0, 7r/2] and y = 0. In either 
case we get matrices of the form 



Aia) = a 1 1 aG [0,1], (5.20) 
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and the matrices that can be obtained from them by interchanging their 
rows and columns. The adjugate of A{a) is 















1-a' 


-(l-a2 





-a -a') 


I -a' 



A'{a)= 1-a^ -(l-a^) 1 (5.21) 

\ -(l-a^) l-a^ J 

The next lemma is an analogue of Lemma 14.11 when the off diagonal 
elements of the kernel are all greater than or equal to zero. It also shows 
that the necessary condition in Lemma [2.61 is satisfied in this case. 

Lemma 5.2 Let A^ be the matrix 

f 1 «'i 4 \ 

A+ = a'2 1 h\ . (5.22) 

V 4 h', 1 ; 

in which a'^, a^ih'i^h^-, d^^ d2 are all greater than or equal to zero, with 
a'^ 02, &1&2, c'^Cg all less than or equal to 1. If A^ is not diagonally equiv- 
alent to a symmetric matrix and A^^ is an M-matrix, there exists a 
diagonal matrix $, with strictly positive diagonal entries, such that $y4+ 
has only one real eigenvalue. 

Furthermore, the real part of the complex eigenvalues of $A+ is less 
than the real eigenvalue. 

Proof of Lemma 15.21 Assume, to begin, that a'^, a'g, 6'^, 62, c'^^, Cg are all 
strictly positive. We use the notation of the proof of Lemma 13.11 except 
that we consider the associated symmetric matrix Kj^ be the matrix 



K+= \ a I h \ , (5.23) 




where a'^a'2 = a^, &162 = b"^ and c'^C2 = c^. It follows from Remark 12.21 
that det i^+ > 0. As in the proof of Lemma 13.11 

det {^K+ - XI) (5.24) 

= -A^ + A^ (01 + 02 + 03) 

-A (0102(1 - a^) + 0203(1 - 6') + 0103(1 - C^)) + 010203 dct K 

= (A -1)2(010203 det K+- A). 
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Note that 



det {^K+ - A/) = -6A + 2(01 + 02 + 03)- (5.25) 



t/A2 
Using the fact that 

(j)i + (j)2 + h > ^ (5.26) 

we see that det ($i^+ — A/) has a local minimum of at A = 1. 
Consider 

det {^A^ - A/) (5.27) 

By Lemma 12.21 it suffices to take A^ to be 







i a c \ 








a I bi 


(5.28) 






V c 62 1 / 




where 6162 '- 


= 6^. We have 








det K+ = 


--1-a^ -b"^ -c^ + 2abc 


(5.29) 


and 










det A+ = l- 


-a^ ~b'^ -c^ + ac(bi + 62). 


(5.30) 



Since 6162 = b"^ we see from the first equality in (15.241) that 

det ($A+ - A/) - det {^K+ - XI) = ac{bi + 62 - 26). (5.31) 

Since det ($-ft'+ — A/) has a local minimum of at A = 1, when 61 7^ 62, 
det ($A+ — A/) has a local minimum that is strictly positive at A = 1. 
This implies that when 61 7^ 62, det ($A+ — A/) has only one real root. 
This is equivalent to the statement of this lemma. 

Suppose now that one of the terms a[, a^, b'l, 62, c[, C2 is equal to zero. 
Let's say it is a[. Then since 



^;' = ^^^ I b[c[ - a'2 1 - c[c', a'A -b[ I , (5.32) 

we see that either 62 or C2 must be zero. Suppose it is 62 . Then A+ in 
(I5.28P has the form 

'1 04 
A+ =14 1 b[ I . (5.33) 

c[ 1 
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and therefore it is equivalent to the symmetric matrix 

A+ = 1 . (5.34) 




We now show that when $A+ has only one real eigenvalue, the real 
part of the complex eigenvalues is less than the real eigenvalue. Let 5 
denote the real eigenvalue. Referring to fl4.17p we see that we must show 
that 

(pi + (p2 + h< 35. (5.35) 

It follows from (I5.3ip that the roots of det ($A_,_ — XI) are greater than 
the roots of det ($K+ — A/). Therefore, by (13. 3p . 6 > 010203 det ii'+. 
Thus to obtain (I5.35P it suffices to show that 

01 + 02 + 03 < 3010203 dct K+. (5.36) 

By (I5.25P the second derivative of \^K^ — XI\ is negative when A > (0i + 
</'2 + 03) /3. Considering the graph of \^K+ — A/| we see that this implies 
that the single real root of \^K^ — A/| is greater than (0i + 02 + 03)/3; 
hence we have (I5.35p . n 



6 Premanenetal vectors with pairwise in- 
dependent components. 

Proof of Theorem 11.21 We first consider the case in which 6 G R^. 
Let r be the kernel of ^. It is enough to consider the case in which 
all the diagonal elements of F are equal to one. Since 6 has pairwise 
independent components we know that (11.61) holds. In this case either 
det r = r(l, l)r(2, 2)r(3, 3) = l or else 




(6.1) 



with abc 7^ 0, or F^ is equal to this matrix. (It is obvious F must contain 
three zeros off the diagonal. Any configuration other than (16. ip or its 
transpose has determinant equal to 1.) 
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Suppose the off diagonal elements of F are positive. It is obvious from 
fl5.2p that r~^, if it exists, is not an M- matrix. Therefore, by Theorem 
l.H r is diagonally equivalent to a symmetric matrix, which implies that 
ahc = 0. If the off diagonal elements of F are negative Theorem 1 1 . 1 1 again 
implies that F is diagonally equivalent to a symmetric matrix, which 
again implies that abc = 0. Therefore Theorem 11.21 holds for 6 G R^. 

We next consider a generalization of (16. ip to nxn matrices, F„, n > 3. 
These are matrices for which 



|/ + aF„| = J](l + a,)+C(F„)J] 



ca 



(6.2) 



i=l 



i=l 



in which C(F„) is a real valued function of the components of F„. We 
have the following lemma: 

Lemma 6.1 Let 9 G R'\, n > 3, he a P -permanental vector, with pair- 
wise independent components and kernel F„. Then C(F„) = 0, in which 
case the components of 9 are independent. 

Proof We show in the beginning of this section that this lemma is true 
when n = 3. For n > 3 the Laplace transform of 9 is 



$(ai. 



ttn) 



n(i+«o+c(F„)n 



on 



i=l 



1=1 



(6.3) 



Let ai = 1 for all 4 < i < n. By Lemma |2T 
$(ai,a2,a3, 1 1 



$(«!, 02,03) 



$(0,0,0,1,...,!) 



;i + ai)(l + a2)(l + a3) + 



(6.4) 



2n-3 



■aia2C(3 



is a Laplace transform of a random variable in R^. Furthermore, the 
form of the right-hand side of (16. 4 p shows that (^1, 92, 9^) is a permanental 
vector with a kernel of the form of (I6.ip . (or its transpose). We show 
in the beginning of this section that for such a vector we must have 
C(F„) = 0. D 
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Proof of Theorem 11.21 continued The proof is by induction. We show 
in the beginning of this section that Theorem 11.21 holds for for n = 3. 
Let n > 3 and assume that the theorem holds for all m < n. Let 6 be 
a /3-permanental process in i?" with pairwise independent components 
and kernel G. This theorem follows from Lemma 16.11 once we show that 
the determinant |/ + aG\ for this process has the form of (16. 2p . 

Suppose it does not. Then the it must contain a term of the form 

k 

C'{G)l[{l + a,^G{i„z,)) (6.5) 

i=i 

where 1 < k < n and (ii, . . . ,ik) is a proper subset of (1, . . . ,n) and 
C'{G) 7^ 0. In fact we know what C'{G) is. Let (a^j, . . . , aj^_j.) be 
the elements of (ai, . . . , a„) that are not in (oj^, . . . , ctj^,) and G' be the 
{n — k) X {n — k) matrix obtained by removing the ii-th,. . . ji^-th row 
and column from G. It is easy to see that 

C\G) = \I + -iG'\ (6.6) 

where / is the {n — k) x {n — k) identity matrix and 7 is the diagonal 
matrix with (7j^j = aj.), 1 < i < n — k. 

By the hypothesis of this theorem \I+'~fG'\~'^ is the Laplace transform 
of a permanental process in K^' with pairwise independent components. 
By the induction hypotheses these components are independent. There- 
fore 

n—k 

C'iG) = l[{l + a,MJ^,J^))■ (6.7) 

1=1 

This shows that the determinant |/ + aG\ for this process has the form 

of dSSD- □ 



7 Proof of Corollary 11.1 



It is easy to see what we must show. For all x,y,z E T 

d{x , y) < d{x , z) + d{y , z) . (7.1) 



28 



This follows if we can show that the determinant of 

/ T{x,x) {T{x,y)T{y,x)y/^ {T{x, z)T{z,x)y/^ 

r= (r(x,i/)r(y,x))V2 T{y,y) {T{y,z)T{z,y)y/^ 

\iTix,z)Tiz,x)y/' iTiy,z)Tiz,y)y/' Tiz,z) 

(7.2) 
is greater than or equal to zero. This follows because, by (12. 2p . the 2x2 
principal minor of F is greater than or equal to zero. Therefore, if |r| > 0, 
r is positive definite and hence the covariance of a Gaussian vector in 
R^. Considering (ITTT]) we get (171]) . 

By Lemma l2^ the kernel of the permanental vector {P{x), P{y), P^z)) 

is 

/ T{x,x) T{x,y) T{x,z) \ 

r= T{y,x) T{y,y) T{y,z) . (7.3) 

Vr(^,x) T{z,y) T{z,z) J 

Since F is the kernel of a permanental vector we have |F| > 0. We 
must show that 

|F| > implies that |f| > 0. (7.4) 

One of the idiosyncrasies that we must take into account is that the off 
diagonal elements of F are always greater than or equal to zero whereas 
the off diagonal elements of F may be negative. 

Suppose F > and |F| > 0. By Theorem 11.11 and Lemma [2.31 either 
F~^ is an M-matrix or else F is diagonally equivalent to a symmetric 
positive definite matrix. In the first case, by Lemma [2741 and Remark l2.2l 
|F| > 0. In the second case, it is easy to see from (12. 151) that |F| = |F| so 
we also have |F| > 0. 

When F has negative off diagonal terms we consider two cases. The 
first is that F is diagonally equivalent to a matrix with positive off diag- 
onal terms, say F', i.e., F = DT'D^^. In this case |F| = \Dr'D^^\ and 
the matrix in (17. 2p is the same for F and F'. Therefore, the argument in 
the previous paragraph shows that |F| > 0. 

Now, suppose F has negative off diagonal terms and it is not diag- 
onally equivalent to a matrix with positive off diagonal terms. Relabel 
the matrix F'. Without loss of generality we can assume 

-T{x,y) -T{x,z) 
r'= I -r{y,x) r{y,y) -r{y,z) I . (7.5) 

-r(2,y) T{z,z) 
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By Lemma 13.11 T' is diagonally equivalent to a symmetric matrix. By 
Lemma 1221 |r'| is equal to 

T{x,x) -{T{x,y)T{y,x)y/' -{T{x,z)T{z,x)y/^ 

-{r{x, y)r{y, x)Y/^ r{y, y) -(r(y, z)T{z, y)Y'^ 

-(r(x,z)r(^,x))V2 -[T{y,z)T{z,y)fl^ T{z,z) 

(7.6) 
Therefore, since F' is the determinant of a permanental vector |r'| > 0. 
This implies that |r| > in this case also. n 



8 M-matrices and infinite divisibilty 

Critical in the proof of Theorem 11.11 is the fact that if a matrix F is 
invertible and F~^ is diagonally equivalent to an M-matrix, then F is the 
kernel of a permanental vector for all /3 > 0. This implies that the vector 
is infinitely divisible. Necessary and sufficient conditions for a vector of 
Gaussian squares to be infinitely divisible are due to R. Bapat and R. 
C Griffiths [U H]. An exposition of this material is given in P^ Chapter 
13]. Eisenbaum and Kaspi, [3^, Lemma 4.2] recognize that the sufficient 
condition in the Bapat-Griffiths criteria for infinite divisibility in the case 
of symmetric kernels also works for non-symmetric kernels. Their proof 
involves probabilistic considerations. Since Theorem 11.11 is only for 3 x 3 
matrices it seems appropriate to give a proof for finite matrices involving 
only linear algebra. 

As one might expect the proof mimics the proof in the symmetric 
case which is given in [H Lemma 14.9.4 and Theorem 13.2.1]. In fact the 
more subtle revisions need to be made in the lemma. We do this next. 
What we prove is given in p:, Theorem 2.3]. The proof in [2] omits all 
details. 

Lemma 8.1 Let B = {Bij}i<ij<p be an invertible matrix such that 
Bi^i > for all i and Bij < for all i ^ j ■ The following are equivalent: 

1. B is an M-matrix; i.e., B^^ > 0. 

2. There exists a positive diagonal matrix D such that all the eigen- 
values of BD have strictly positive real parts. 
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3. There exists a matrix C > such that BD + C = XI where A is 
greater than the spectral radius of C. (I.e., \ > p, where p := p{C) 
is the magnitude of the maximal eigenvalue of C .) 

Proof We first show that 1. imphes 2. Let D be a diagonal matrix 
with diagonal entries Di^i = J2^=i{^~^}i,j^ i = 1, . . . ,n. Since B~^ > 0, 
BD has positive diagonal elements and negative off-diagonal elements. 
Note that 

2_^{BD}ik = 2_^^i,kDk,k = 2_^^hk2_^{^}k,j (8-1) 

k k k j 

= ^ Bi^k{B~^}k,j = ^ Sij = 1. 

This shows that that BD is strictly diagonally dominant, i.e., for all i, 

J2 |5A,fc| < 5A,.. (8.2) 

k:k^i 

Consider the symmetric matrix A = BD + {BDy, where the super- 
script t represent conjugate transpose. To see that A is positive definite, 
let u = (mi, . . . ,Up) be any real vector, and let Qj := {BD}ij. Using 
([H2D we have 

{BDu,u) = y^u'^Ci^i+ y^ UiUjCij (8.3) 

It follows from this that 

{Au,u) > ^ [uf + u] ~ 2\ui\\uj\) (IqjI + \cj^i\) 



i,j:i<j 



14) 
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Since A is a symmetric and strictly positive definite matrix its eigenvalues 
are all greater than 0. 

In general, if a is an eigenvalue of BD^ a is an eigenvalue of {BDy. 
Let A be the eigenvalue matrix of BD, and write BD = EAE^, where E 
is the matrix of eigenvectors of BD. Then 

A = EAE^ + EA^E^ = E{A + A^)E\ (8.5) 

Therefore, since the eigenvalues of A are all greater than 0, the eigenval- 
ues of BD have strictly positive real parts. 

To show that 2. implies 3. we first note that because the off diagonal 
elements of BD are less than or equal to zero, it is easy to see that we 
can find a positive number A and matrix C" > such that 

BD = XI - C. (8.6) 

It remains to show that we can do this with A > p{C'). 

Clearly, for any p-dimensional complex valued unit vector y, \C'y\ < 
p. Let X be an eigenvector of C", of unit length, such that C'x = •yx, 
where I7I = p. Let x be the vector with Re{xi) = \Re{xi)\ and Im{xi) = 
Im{xi) for all i. Trivially, |a;| = \x\. Then, since C is positive, un- 
less X = X, \C'x\ > \C'x\ = p. Here we have assumed that a least 
one real component of x is not equal to 0. If that is not the case then 
C'{ix) = '~f{ix) and since ix has real components we can apply the argu- 
ment above to see that we must have x = ix. Consequently, there exists 
an eigenvector x of C such that Re{xi) > for all 1 < i < p, with strict 
inequality for at least one i, satisfying 

C'x = 7X. (8.7) 

Let A be the strictly positive definite symmetric matrix defined above 
in this proof. The argument used to obtain (18. 5p shows that there exists 
a A such that 

2A>7 + 7. (8.8) 

Let 7 = a + z6, in which a and b are real. By (18. 8p there exists an e > 
such that X = a + e. 

For any 5 > 0, BD + {C + 61) = (A + 6)1. The spectral radius of 
C + 61 is \{a + 6)"^ + fe^l-*^/^. It is easy to see that for sufficiently large 6, 
say 60, this is smaller than (A + 60)"^. Taking C = {C + 6qI) we get 3. 
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We conclude the proof by showing that 3. imphes 1. Clearly the 
absolute values of all the eigenvalues of C /\ are less than one. Therefore 
we can write 

to see that (BD)^^ = D^^B^^ > 0. Since D is a positive diagonal 
matrix, 1. follows by multiplying D^^B^^ by D. n 

In the next theorem, rather than focus on infinite divisibility, we 
characterize the kernels that give /3-permanental processes for all /3 > 0. 

Theorem 8.1 Let B be an n x n invertible matrix and let A = B^^ : 
If there exists a signature matrix M such that MBM is an M-matrix, 
then Man is a kernel of a (5-permanental vector for all (5 > 0. (I.e., 
\I + aA/'AA^|~^ is the Laplace transform of a random positive vector for 
all 13 >0.; 

Proof One of the main steps in this proof is that for any matrix U with 
spectral radius less than 1, 

iog(i/-tn)=f:^i5£!E!}, (8,10, 

^-^ n 

n=l 

This is explained on [H| page 562], for matrices that are closely related to 
symmetric matrices. However the proof goes through exactly as written 
for any matrix with the spectral radius less than 1. 

We proceed with the proof. To simplify the notation we might as well 
assume that B itself is a M-matrix, i.e. that we can take M = I . Then 
if we can show that there exist two strictly positive diagonal matrices D'^ 
and D'2 such that D'^AD'2 is a kernel of a /3-permanental vector for all 
/3 > 0, the conclusion of 1. follows from Lemma [2. 1[ 

Since B is an M-matrix, it follows from Lemma [8. ![ 3. that for some 
positive diagonal matrix D, B = DB = XI ~ C for some positive matrix 
C and A with the property that A > p{B). Let A = {B)~^ and 5* be a 
diagonal matrix with diagonal entries 0<Sj<l,i = l,...,n. Then for 
any diagonal matrix a with diagonal entries a^, i = 1, . . . , n, we can can 
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find a positive number t and a diagonal matrix 5* such that a 
We have 



\I + aA\ 



\I + t{I-S)A\ 
\A\\B + t{I-S) I 
\A\\{X + t)I-iC + tS)\ 



tiI~S). 



(8.11) 



(A + t) 

Since S has all entries less than 1 and all eigenvalues of C are less 
than A in magnitude, it follows from (18.101) that 



log(|/ + aA\) = log{A) + log(t + A) + ^ 



fc=i 



trace{(C + t5)^"} 
k{XTW^ 



(8.12) 



for sufficiently large A > 0. 

Since C is positive it is obvious that the expansion of this power series 
in about S = has non-negative coefficients. By [H Lemma 13.2.2], this 
implies that |/ + a^l"'^ is the Laplace transform of a positive random 
vector for all /3. Since A = D~^B, the comments in the second paragraph 
completes the proof. n 



Remark 8.1 ^0\ Proposition 4.6] states that a necessary condition for 
a kernel F to define a permanental vector in (II. ip . for all /3 > 0, is that F 
and all matrices obtained from F by multiplying its rows by non-negative 
real numbers numbers, have only real non-negative eigenvalues. This is 
not correct even for symmetric matrices since kernels of Gaussian squares 
that are in class 1. and not in class 2. satisfy this condition but, since 
they are not in class 2., they are not kernels of permanental vectors for all 
P > 0. In fact, considering Lemmas 14.11 and 15. 2[ a correct statement is: 
If a 3 X 3 matrix F is the kernel of a permanental vector, and it, and all 
matrices obtained from it by multiplying its rows by non-negative real 
numbers have only real non-negative eigenvalues, then F is diagonally 
equivalent to a symmetric positive definite matrix. 
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